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ABSTRACT 

Using the "Scenario Machine" we have carried out a population synthesis of X-ray binaries 
for the purpose of modeUing of X-ray luminosity functions (XLFs) in different types of galaxies: 
star burst, spiral, and elliptical. This computer code allows to calculate, by using Monte Carlo 
simulations, the evolution of a large ensemble of binary systems, with proper accounting for the 
spin evolution of magnetized neutron stars. 

We show that the XLF has no universal type. It depends on the star formation rate in the 
galaxy. Also it is of importance to take into account the evolution of binary systems and life 
times of X-ray stages in theoretical models of such functions. We have calculated cumulative and 
differential XLFs for the galaxy with the constant star formation rate. Also we have calculated 
cumulative luminosity functions for different intervals of time after the star formation burst in 
the galaxy and curves depicting the evolution of the X-ray luminosity after the star formation 
burst in the galaxy. 

Subject headings: binaries: close — binaries: general — X-rays: binaries — X-rays: general 



INTRODUCTION 



The evolution of t he X-ray luminosity o f gala x- 



ies was predicted by iTatarinzeva et ah 
They studied X-rays come from X-ray binary 
stars. The evolution of the X-ray luminosity Ls{t) 
was calculated assuming the (5-function shape for 
the star formation rate (simultaneous birth of the 
stars). The evolution of the X-ray luminosity of 
the galaxy with an arbitrary star formation rate 
(j){t) can be represented as 



ter i > 2 X 10^ years (long time scale) from the 
star formation burst can b e well fitted by power 
law (jTatarinzeva et al. |[l989il ^: 



L{t) « 3-10 



40 



f N \ f t 



-1.56 



erg-s 



(2) 



L{t) = / L&{t-T)<P{T)dT, (1) 

J —CO 

The evolution of the total X-ray luminosity af- 



^See also I Lipunov et al. 



here N is the total number of the stars in the 

ga laxy. 

iLipunov et al. ( 1996a ) studied the evolution of 
stellar populations after the star formation burst 
occurring in the conditions similar to the Milky 
Way, in the central part of the galaxy, on a 
timescale of 10 Myr. Their results include a num- 
ber of X-ray transients (each consisting of a neu- 
tron star and a main sequence star), super ac- 
creting black holes, and binaries consisting of a 
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black hole and a supergiant, as functions of time. 
They showed that absolute and relative numbers 
of massive binary systems including neutron stars 
and black holes can serve as a goo d indicator of 
the ag e of the star formation burst. IPopov et al. 



dN 
dL 



38 



{3.3tl[l)SFR X L~i-''^' '^ L < L„ (4) 



(|1998I ) also made fits to dependencies N{t) for dif- 
ferent types of objects, where N{t) is the number 
of sources, t is the time after the star formation 
burst. 

( 2000[ ) combined 



where L38 — erg s ^ and SFR is measured 

in units of Mq per year. 



Zezas et al. 



Van Bever fc Vanbeveren 



their close binary population number synthesis 
code with the formation mechanism of X-radiation 
in young supernova remnants and in high mass X- 
ray binaries. They demonstrated that the impact 
of interacting binaries is substantial. 

Numerous point- like extragalactic X-ray sources 
were discov ered during last yea r s due to Chan 



(j2004^ presented the X-ray lumi- 
nosity function of the Antennae galaxies based on 
8 observation performed with Chandra, 7 of them 
were obtained between January 2001 and Novem- 
ber 2002. After combining all observations they 
detect a total of 120 sources down to a limiting lu- 
minosity of ~ 2 • 10^^ erg s"-'^. Authors suggested 
that comparison between the XLFs of the individ- 
ual observations showed that they are not affected 
by source variability. The cumulative XLF of the 
coadded observations was represented by a single 



dra (s ee e.g. iMuno et al. I (|2004l ). iGrindlay et al. 



2002 



power law N{> L) ^ L 



-0.52t°\ 



(E005h) and XMM - Newto n (see e.g. iKong \ (l2003h. indication for a 'bump' at ^ lO^^ erg s 



There was an 
but at 



Georgakakis et al. l(|2004 ). lGeorgantopoulos et al. 
j 20051) ) mis si ons. Some auth o rs ( Grimm et al 



20031 : iGilfanov I 12004 iKim k Fabbiano 



20041 ) report about power law X-ray luminosity 
function: 



this point its significance was not clear. If this 
bump is statistically significant it could be evi- 
dence for Eddington limited accretion on compact 
objec ts or anisotropic em ission from the accretion 
disk (jZezas fc Fabbiano | [20G2). 

Belczvnski et al. ( 20041) constructed synthetic 



dN 



X SFR,a^ 1.5, 



(3) 



where SFR is the star formation rate. 

These data were discussed by IPostnov I (|2003l ) 
from theoretical point of view. 

Grimm et ahl (|2003h realized that, within the 
accuracy of the presently available data, a lin- 
ear relation between high mass X-ray binaries 
(HMXB) number and star formation rate (SFR) 
exists. They suggest that the relation between 
SFR and collective luminosity of HMXBs is non- 
linear in the low-SFR regime, ~ SFR^^ ''^ , and 
becomes linear only for a sufficiently high star for- 
mation rate, SFR > 4.5Mq yr'^ (for M > 8Mq). 
Also they obtained the universal luminosity func- 
tion of HMXBs and fitted the combined luminos- 
ity function of M82, Antennae, NGC 4579, 4736 
and Circinus using a maximum-likelihood method 
with a power law with a cut-off at = 2.1 • 10'*° 
erg s~* and normalized the result to the combined 
SFR of the galaxies. Their best-fitting luminosity 
function in the differential form is given by 



X-ray binary populations for direct comparison 
with the X-ray luminosity function of NGC 1569 
observed with Chandra. They produced hybrid 
models meant to represent the two stellar popu- 
lations: one old and metal-poor, with continuous 
star formation for ~ 1.5 Gyr; and another a recent 
and metal-rich population. They found that for 
typical binary evolution parameters, it is possible 
to quite closely match the observed XLF shape. 

Our critical points concerning both observa- 
tional and theoretical aspects is in that that there 
is no observed universal luminosity function be- 
cause: 

1. Number of bright X-ray binaries is very 
small per galaxy. 

2. We do not know real X-ray luminosity due 
to high variability of binary X-ray sources, 
on scales from seconds up to 100 years. 

There is no simple (with one slope) theoretical 
universal luminosity function because: 

1. X-ray population is the mix of different types 
of binaries with different mass exchange 
types. 



2 



2. Number of the systems with definite lumi- 
nosity depends on spin evolution of a neu- 
tron star which has no direct connection to 
mass of its companion. 

3. Theoretical arguments for universal function 
being at present time are not quite correct, 
because they exclude life-times (which de- 
pend on optical compan ion mass] of bi nary 
stars in accretion stage ( Postnov"1l2003f l. 



We stress that it is of great importance to take 
the spin evolution of NSs into account. Quantity 
of accreting neutron stars which give their con- 
tribution to the luminosity function is determined 
by their magnetic fields and spin periods. Neu- 
tron stars can be in a non-ac creting state (pr o- 



peller, ejector, see for details iLipunov I mm). 



This circumstance usually is not taken into ac- 
count in population synthesis models. 

We must observe much more sources and deter- 
mine their types to make correct luminosity func- 
tion. In any case XLFs must have different slope 
for different types, ages and star formation histo- 
ries in galaxies. 

Ultra luminous X-ray sources (ULXs) with 
Lx > lO^^erg s~^ have been discovered in great 
amounts in external gal axies with ROSAT, Chan - 

l|200fih 



dra and XMM-Newton. iRappaport et al 



carried out a theoretical study to test whether 
a large fraction of the ULXs, especially those in 
galaxies with recent star formation activity, can be 
explained with binary systems containing stellar- 
mass BHs. To this end, they have applied a unique 
set of binary evolution models for BH X-ray bi- 
naries, coupled to a binary population synthesis 
code, to model the ULXs observed in external 
galaxies. They find that for donor stars with ini- 
tial masses > IQMq the mass transfer driven by 
the normal nuclear evolution of the donor star is 
sufficient to potentially power most ULXs. This 
is the case during core hydrogen burning and, to 
an even more pronounced degree, while the donor 
star ascends the giant branch, although the latter 
phases last only 5 per cent of the main-sequence 
phase. They show that with only a modest vio- 
lation of the Eddington limit, e.g. a factor of 10, 
both the numbers and properties of the majority 
of the ULXs can be reproduced. One of their con- 
clusions is that if stellar-mass BH binaries account 
for a significant fraction of ULXs in star-forming 



galaxies, then the rate of formation of such sys- 
tems is 3 • 10~^ yr~^ normalized to a core-collapse 
supernova rate of 0.01 yr~^. 



King et al. 



( 20011 ) investigated models for the 
class of ultraluminous non-nuclear X-ray sources 
(ULXs) seen in a number of galaxies and probably 
associated with star-forming regions. The assump- 
tion of mild X-ray beaming suggests instead that 
ULXs may represent a short-lived but extremely 
common stage in the evolution of a wide class of 
X-ray binaries. The best candidate for this is the 
phase of thermal-timescale mass transfer that is 
inevitable in many intermediate- and high-mass 
X-ray binaries. This in turn suggests a link with 
the Galactic microquasars. The short lifetimes of 
high-mass X-ray binaries would explain the asso- 
ciation of ULXs with episodes of star formation. 
These considerations still allow the possibility that 
individual ULXs may contain extremely massive 
black holes. 

We also would like to remember the old con- 
sideration of the supercritical non-sph erical ac 



cretion onto magnetized neutron stars (jLipunov 



1982al lbh. In this case matter falls to the magnetic 
poles of the neutron star. Maximum energy re- 
lease proves to be i = ^^LEddip-'io)'^^'^ , where /iao 
- is the magnetic dipole moment of the neutron 
star in 10^° G cm^. 

2. DESCRIPTION OF MODELS 

2.1. Binaries under consideration and 
mechanisms of mass accretion 

The "Scenario Machine" code includes the next 
types of the mass accretion by the compact star: 

1. Accretion from the stellar wind. 

2. Accretion from the disk-shaped stellar wind 
of Be-stars. 

3. Mass transfer through the inner Lagrange 
point during Roche lobe overfiow stage: 

(a) On thermal timescale. 

(b) On nuclear timescale. 

(c) On magnetic stellar wind timescale. 

(d) On gravitational radiation timescale. 
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Induced stellar wind is not included into the 
program. 

Most of the X-ray pulsars in the Milky Way be- 
long to binaries which consist of the Be-star and 
the neutron star j Liu et al. I bOOOHLii et al. 112001 : 
Raguzova fc Popov 2005). The mass loss by the 



Be-star is determined by its rotation. Its mass 
influences onto its wind to a lesser degree. At 
the same time we see a little part of the X-ray 
binaries consisting of Be- and neutron stars due 
to variability of mass transfer processes and the 
transient character of accretion in such systems 



(jvan den Heuvel 111994 ). 



So, we should study as much types of X-ray 
binaries as possible. This is the list of the systems 
under our consideration: 

1. NA-I-I: the accreting neutron star with the 
main sequence companion. 

2. NA-I-II: the accreting neutron star with the 
super giant companion. 

3. NA+III: the accreting neutron star with the 
companion filling its Roche lobe. 

4. NA-t-Be: the accreting neutron star with the 
Be-star companion. 

5. BH-l-II: the black hole with the super giant 
companion. 

6. BH+III; the black hole with the companion 
filling its Roche lobe. 

7. SNA+III: the super accreting neutron star 
with the companion filling its Roche lobe. 

8. SBH-I-III: the super accreting black hole 
with the companion filling its Roche lobe. 

The last two types of systems are taken into 
consideration for the purpose of modelling of 
ULXs. Radiation o f such objects can be strongl y 



coUimated (see e.g. ICherepashchuk et al. I ( 20051 ) 1 
to a degree ~ 1°. Also we take into account pos- 



sibility of mild beaming (see e.g. iKing et al. 
(|200l[ )). If the radiation of the source is coUi- 
mated, then we should reduce calculated number 
of binaries using formula 



No = —N,, 
An 



because we can not see the object if its beam is 
directed away from us. We recalculate X-ray lu- 
minosity of such systems using formula 



_4n 



(6) 



in order to obtain the luminosity under the formal 
assumption of spherically symmetric radiation. In 
these equations 51 is the doubled solid coUimation 
angle of the radiation, Lc is the calculated lumi- 
nosity of the source and Nc is the calculated num- 
ber of sources, Lo and No are the same observable 
values. 

We have to say some words about Wolf-Rayet 
(WR) stars with black holes or neutron stars. 
Number of binaries consisting of the accreting 
black hole and the WR-star is very small, because 
WR-stars have very high velocity of the wind. 
That is why accretion disks are not able to form 
in wide pairs (with orbital period > 10 hours; or- 
bital period of Cyg X-3, for e xample, is ~ 5 hour s; 



see fo r detailed description iKarpov fc Lipunov 
There are no binaries consisting of WR- 



stars and accreting NSs, because NSs accelerate 
their rotation during second mass exchange (recy- 
cling) and theref ore become propellers or ejectors 
(jLipunov Ill982a l in such kind of binaries. 

Note that our conclusions concerning accreting 
neutron stars with Be-stars, super accreting neu- 
tron stars with non-degenerate stars filling their 
Roche lobes, super accreting black holes with non- 
degenerate stars filling their Roche lobes have ap- 
proximate character, because it is impossible to 
depict correctly temporal and angular dependen- 
cies of their radiation at present time. Our calcu- 
lations show that real luminosity function is com- 
pound. 

2.2. List of main evolutionary parameters 

Since the algorithms used in the "Scenario Ma- 
chine" have been described many times, we shall 
only note the most important evolutionary pa- 
rameters influencing the results of the numeri- 
cal modeling of the binaries under consideration. 
A detailed description of the "Sce nario Machine" 



can be f ound i n the next works: iLipunov et al. 
(Il996bi r 



20071 ). 



(5) 



The initial masses of primary components were 
varied between IOM0 and 120Mo. To describe 
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also a kind of ULX objects consisting of a black 
hole with mass higher than ~ IOOMq and an op- 
tical star in Roche lobe overflow stage we have 
conducted a population synthesis also with lower 
and upper limits equal to 120Mq and lOOOM© 
correspondingly. 

We assume zero initial eccentricity, all initial 
mass ratios have equal probability, initial mass of 
the secondary star is in the range O.IMq - mass 
of the primary star. 

Mass loss by optical stars in the course of 
their evolution remains incompletely explored. 
Despite the fact that is has been possible to 
substantially reduc e the uncertainties (see, e.g., 



(|Bogomazov et al. 



2005[) ). no clear justifica- 
tion for a choice of a standard scenario has 
emerged. Therefore, we carried out our compu- 
tations for two scenarios for mass loss by non- 
degenerate stars, which we call A and C. A de- 
tai led description of the se models can be found 



in (jLipunov et al. 1120071 ). Scenario A has a weak 



stellar wind. The stellar wind of massive stars 
(with masses higher than 15Mq) is higher in sce- 
nario C, for lower-mass stars, scenarios A and C 
are equivalent. The total mass loss in any evo- 
lutionary stage is higher in scenario C than in 
scenario A. 

Common envelope stage efficiency ac e is equal 
to 0.5. 

Minimal initial mass of the star which produces 
a black hole as the result of its evolution is 25Mq. 
We assume the par ameter khh = Mhh I M Pre.fi n to 



(|200a for 



be equal to 0.5 (see iBogomaz ov et al. 
detailes), MpreSN is the mass of the pre-supernova 
star which produced the black hole with mass M^h- 

Initial mass of the new-born neutron star is ran- 
domly distributed in the range 1.25 - IAAMq. 
Maximum mass of the NS (Oppenheimer-Volkoff 
limit) equals to Mqv = 2.0Mo in our calcula- 
tions. Initial value of the magnetic field of NSs is 
assumed to be equal to 10^^ Gs, the field decay 
time is assumed to be equal to 10^ years. Charac- 
teristic kick velocity of the new-born neutron star 
we accept to be equal to 80 km s^^ in this work. 

We use two different values of collimation angle 
for supercritical regimes of accretion: a = 1° and 
a = 10°. 



2.3. Result normalization 

Birth frequencies of binaries were calculated us- 
ing the next formula: 



1 



(7) 



here Vgai is the frequency of birth of a specific bi- 
nary system type in a spiral galaxy, Ncaic is the 
number of the systems under our consideration ap- 
peared during calculations, Ntr is the total num- 
ber of binaries which evolution was calculated, Mi 
is the minimal initial mass of a star in our calcu- 
lations. We treat a spiral galaxy in this case as a 
galaxy with constant star formation rate which is 
defined by the Salpeter function. 

Quantities of the systems in a spiral galaxy were 
calculated using equation ([5]). 



Ntr Ml 



1.35 ' 



(8) 



here Ngai is the quantity of a specific binary sys- 
tem type in a spiral galaxy, ti is the life time of 
the binary system under consideration. 

Along with modeling population in the spiral 
galaxy we also made some estimations of evolu- 
tion of X-ray luminosity function and total X-ray 
luminosity in an elliptical galaxy. Quantities of 
the systems in the elliptical galaxy were calculated 
using equation 



Nnnl = K, 



M, 



gal 



ScM 



Ml 



ScM 



Igal 



-1.35 



AT ' 



(9) 



here Mgai = lO^M© is the mass of typical galaxy, 
MiscM is the minimal initial mass of a star in 
our calculations, Migai — O.IMq is the minimal 
initial mass of a star, ti is the life time of a binary 
system under our consideration in the range of the 
time interval between T and T-f- AT. We treat an 
elliptical galaxy in this work as a kind of object 
in which all stars born at the same time and then 
evolve ((5-function star formation rate). 

2.4. Constraints on key parameters of the 
evolutionary scenario 

Previous estimates of the ranges of parame- 
ters deter mining the evolution of binari es were ob- 
tained by iLipunov etaTl (|l996d.ll997l ). 
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Since that time, some new results related to the 
evolution of binaries have been obtained, and we 
carried out additional computations of constraints 
that can be applied to the parameters of the evo- 
lutionary scenario. The latest observational es- 
timates of the kick velocities receiv ed by NSs in 
super n ovae explosions are given by iHobbs et al. 



(|2005l ). where it is concluded that the typical kick 
magnitude is cr — 265 km s^^. 

An attempt to obtain a more accurate estimate 
of the mass- loss effici ency in the comm on-envelope 
stage was made by IPewi k. Tauris I (|2000.) . who 
tried to take into account the concentration of the 
stellar material toward the center: j^^^"^"'' • 
However, they assumed that the efficiency in the 
common-envelope stage was ^ce = 1. In general, 
this parameter is not known accurately. Our coef- 
ficient ac E is the p roduct of Mr? f, and the p arame- 



ter A considered bv lPewi fc Tauris ( 200Clf ). which 



describes the concentration of the stellar matter 
toward the center. For this reason, we use the 



value of ac_B suggested bv lLipunov et al. I (|l996cl ). 

We would like to note one more important cir- 
cumstance. Ill-defined parameters of the evolu- 
tionary scenario, such as fo, oicE, the stellar wind 
efficiency, and so on, are internal parameters of the 
population synthesis. In the future, they may be 
defined more precisely, or their physical meaning 
may change: the kick-velocity distribution may 
turn out not to be Maxwellian, it may be that 
the complex hydrodynamics of common envelopes 
cannot be described using the parameters ace and 
A, the mass ratio distribution f{q) may be not a 
power law. There exists only one way to verify our 
results: comparison of our model predictions with 
observational data. 

For this reason, we suggest two quantities to be 
compared to test the model: the ratio of the calcu- 
lated and observed numbers of Cyg X-3 systems, 
and the ratio of the number of binary radio pul- 
sars with NS companions and the total number of 
radio pulsars (both single and binary), ^^^'''"'^ , 
where Npsr+NS is the number of radio pulsar in 
binary systems with a neutron star, Npsr is the to- 
tal number of radio pulsars, binary and single. To 
avoid the need to differentiate between young pul- 
sars and old pulsars that have been accelerated by 
accretion, we consider only young pulsars. Note 
that the observed value of this ratio is ^ 0.001 
( ATNF catalogue I l200(ih : among more than 1500 



known single and binary radio pulsars, only two 
young pulsars in pairs with NSs have been discov- 
ered (J2305-I-4707 and J0737-3039). As a model 
of a Cyg X-3 type system, we take a binary con- 
taining a BH with WR companion that is more 
massive than > IMq and having an orbital pe- 
riod 10 hours. 

Figure [3] of the calculated number of binaries 
with a NS and radio pulsar Npsr+NS and the cal- 
culated sum of the numbers of single and binary 
radio pulsars Npar depends on the kick velocity 
wq. The width of the shaded region reflects the 
variation of the efficiency of the common envelope 
stage acE in the range 0.21.0. 



Figure[3]shows the OCCO criterion (jLipunov et al. 
Il996bl ) for the ratio ^^°''+"^' . The typical kick ve- 



locity vq is plotted along the horizontal axis. The 
width of the shaded region reflects the variation 
of the efficiency of the common envelope stage 
acE in the range 0.21.0. The observed value of 



is ~ 0.001. 



As seen from Figs. [3] and [21 the characteristic 
value of kick velocity vq cannot exceed « 200 km 
s~^. B y this reason we mak e use of the results of 
paper (jLipunov et al. |[l997l ). 

Figure [3] shows the number of Galactic Cyg X-3 
systems in our model as a function of the common 
envelope efficiency. This figure shows that we can 
essentially exclude values acE < 0.3. 

3. RESULTS AND CONCLUSIONS 

Four simulations runs were performed, each 
simulating the evolution of 1 ■ 10^ binary sys- 
tems. Two of them were performed with weak 
stellar wind (stellar wind type A), and other mod- 
els with reasonably high stellar wind (stellar wind 
type C). In each of these cases we made our cal- 
culations using two value areas of initial mass of 
the primary star in Salpeter's power law: in the 
range 10 — I2OM0 for investigations of all types 
of systems under consideration, and in the range 
120 — IOOOM0 to qualitatively depict only ultra lu- 
minous objects consisting of super accreting inter- 
mediate mass black holes with companions filling 
their Roche lobes. 

In the Figures [3] and [5] we show birth frequency 
of different types of X-ray sources in the spiral 
galaxy. In the Figures [3] and [7] we present cu- 
mulative luminosity functions of different types of 
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X-ray sources in the same galaxy. Figures [3] and 
[3] were calculated using stellar wind type A (weak 
stellar wind) . Figures [5] and [7] were calculated us- 
ing stellar wind type C (moderate stellar wind). 
Marks in Figures 4-7 are (see abbreviation in 
Section [231) : 1, NA-fl; 2, NA-hll; 3, NA-hlll; 4, 
NA-fBe; 5, BR+ll; 6, BH-HlII; 7a, SNA-I-III, col- 
limation angle (for super critical regimes of accre- 
tion) a ^ 10°; 7b, SNA-^III, a = 1°; 8a, SBH+III, 
a = 10°; 8b, SBH-hlll, a = 1°; 9a, SBH-I-III, 
a = 10°; 9b, a = 1°. For curves 9a, 9b minimal 
initial mass of the primary star is I2OM0, in other 
cases it is equal to IOMq. 

As one can see from Figures [3] - [71 different 
types of X-ray binary systems belong to different 
luminosity ranges, their luminosity functions have 
different slope. These facts are evidence of com- 
plexity of the X-ray luminosity function. 

Comparisons between figures 4 and 5, 6 and 
7 convince us of the importance of taking into 
account life times of X-ray stages in theoretical 
models of XLFs. Relative abundances of differ- 
ent types of X-ray binary systems in the birth 
frequency function and in the luminosity function 
are different. For example, we can see from Fig- 
ure 4 that the birth frequency of NA-I-II X-ray 
binaries is about ten times higher than the birth 
frequency of NA+I X-ray binaries. But the super 
giant life time is much shorter than the life time of 
the main sequence star, so, as we can see from Fig- 
ure 5, quantity of NA+I binaries is only two times 
less than quantity of NA-I-II systems in the spiral 
galaxy. Stronger stellar wind (type C) makes this 
difference even greater (compare Figures 6 and 7) . 

The stellar wind magnitude essentially influ- 
ences the scenario for two reasons. First, the 
spherically symmetric wind leads to increase in 
component separation. Secondly, stellar wind 
greatly affects the final evolutionary outcome of 
massive stars. In particular, the choice of wind 
strength will change the mass d istribution of black 



holes seen in the population (jBogomazov et al. 



20051 ). as the black hole progenitor loses a differ- 
ent amount of mass prior to collapse. Moreover, 
the total mass loss of a star by wind may cause a 
change in its remnant type (it may produce a neu- 
tron star instead of a black hole). We can see from 
Figures 4-7 that stronger stellar wind (type C) 
dramatically decreases quantities of many types of 
X-ray binaries (and affects all kind of them) . 



In the Figures [8] and [9| we show cumulative lu- 
minosity functions of all investigated systems in 
the spiral galaxy like the Milky Way. See Tables 
1 and 2 for numerical data. In these Figures a 
is the collimation angle in supercritical regimes of 
accretion. Figure [8| was calculated under the as- 
sumption of stellar wind type A, Figure[9|was cal- 
culated under the assumption of stellar wind type 
C. 

Figures [5| and [5| show that the X-ray luminos- 
ity function has different slope in different ranges 
of luminosi ty (see also Ta b les 1 and 2 for numer- 
ical data). iGrimm et al. (2002) argued that the 
best values of the slope and normalization of the 
cumulative form of the luminosity function is 



iV(> L) = 5AxSFR(L 



-0.61±0.12 



210 



-0.61±0.12 



) 

(10) 

over the luminosity range between ^ 10'^^ erg s~^ 
and ~ 10^° erg s~^ (see Figure 5 and Equation 7 
in their paper, but they gave narrower luminosity 
range as the result in Conclusions of the article). 
Our calculations show similar XLF slope over the 
ranges between w 2 • 10'^^ erg s"^ and w 10'^* erg 
s~^, and between « 2-10'^^ erg s~^ and « 10'*^ erg 
s^^ (the last range depends on our assumptions 
about the collimation angle of the X-ray emission 
in cases of super critical accretion) . Between these 
two ranges the XLFs (our theoretical curves) be- 
come very steep due to Eddington limit (there are 
a lot of NA-I-III systems, and their Imninosity is 
about this value, see Figures [3| - [7|) . 

In the Figures [TUj and [11] we show birth fre- 
quency (a) of all investigated systems (differen- 
tial function) and differential luminosity function 
(b) of X-ray binary sources in the Galaxy. Marks 
in the Figures are: 1, collimation angle (for su- 
per critical regimes) a = 10°; 2, a = 1°. Figure 
[TOl was calculated under the assumption of stel- 
lar wind type A, Figure [Tl] was calculated under 
the assumption of stellar wind type C. Luminos- 
ity functions in differential form also have different 
slope, there is no evidence for the universal XLF. 

In the Figure[12|we show cumulative luminosity 
functions of all investigated systems in the ellip- 
tical galaxy after the star formation burst. The 
curves in the Figure represent the next models: 
1, stellar wind type A, collimation angle (for su- 
per critical regimes) a = 10°; 2, wind A, a — 1°; 
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3, wind C, a = 10°; 4, wind C, a = 1°. The 
time ranges after the star formation burst in the 
Figure are: a, 0-10 miUion years; b, 10-100 mil- 
Uon years; c, 100 miUion - 1 biUion years; d, 1-10 
bilHon years, v Figure [T^ shows the evolution of 
luminosity function on a long timescale after the 
stellar formation burst in the elliptical ^ galaxy. 
As one can see from this Figure, there is no ev- 
idence for the universal XLF. Nevertheless, note 
that numbers of systems in this figure are quite 
relative. Any systems were added to the number 
of the systems in appropriate interval of time if 
they born as X-ray system or still show itself as 
X-ray source during this period or part of it, but 
life time of a system can be less than the dura- 
tion of the period and a system can born not in 
the beginning of the period. For more precision 
it is necessary to take less intervals of time, but 
our purpose is to show long time evolution quali- 
tatively. 

iBelczvnski et al.~ (2004) found that the depen- 
dence of the XLF slope on age is non-monotonic in 
the dwarf (post)starburst galaxy NGC 1569. They 
studied behavior with time of theoretical normal- 
ized XLFs for two stellar populations: one old at 
1.5 Gyr and one young at age 10, 70, and 170 
Myr (continuous SFR through 1.5 Gyr, and 10, 
70, and 100 Myr, respectively). The average SFR 
in the old population was assumed to be 20 times 
smaller than that in the young population. Di- 
rect comparison between our results is difficult, 
because we use different star formation models in 
our calculations. One of the common features is in 
that that the XLF should evolve with time. Also 
we suggest that their XLFs can be fitted by bro- 
ken power laws, Belcz ynski et al. 1 (2004) did not 
obtain uniform XLF in NGC 1569. 

In the Figure [13] we show the evolution of the 
X-ray luminosity after the star formation burst 
(T = 0) in the galaxy with mass lO^M©. See 
Table 2 for numerical data. In this Figure: 1, 
our calculations, s tellar wind type A ; 2, th e re- 



sult obtained by iTatarinzeva et a~ I (|l989l) : 3, 



our calculations, stellar wind type C. We should 
note that in comparison w i th res ults obtained by 
Van Bever fc Vanbeveren ( 2000f l we do not take 



into account the X-ray emission from supernova 



remnants in our models. Our data in this Figure 
start at their end point (10 Myr). After 4-10^ mil- 
lion years since star formation burst in the galaxy 
its X-ray luminosity can be rather well fitted by 
power law L{T) ~ T~°; a is equal to 1.56 and 1.8 
in very wide ran ge of time (see Table 3 for details). 



Previous work (jTatarinzeva et aL 19891 ) showed 



approximately the same result which we can con- 
firm. The cause of differences is in that that 16 
years ago calculations were conducted if authors 
were taking into consideration not so much types 
of systems as in the present work. Also models of 
evolution of binaries have changed. Stronger stel- 
lar wind (see Tab le 4) makes the our result almost 
inconsistent with Tatarinzeva et al. ] §989). 
So, our calculations show the next results: 

1. X-ray luminosity function of binary X-ray 
sources is complicated, it has different slope 
in different ranges of luminosity. So, there 
is no universal X-ray luminosity function of 
binary X-ray sources. 

2. X-ray luminosity function of binary X- 
ray sources depends on the star forma- 
tion rate as it was firs t shown in 1989 
( Tatarinzeva et al. 1 ll989h . 



3. It is necessarily to take into account spin 
evolution of neutron stars and life times of all 
stages during theoretical modelling of X-ray 
luminosity function of binary X-ray sources. 
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Fig. 1. — Calculated number of Cyg X-3 type sys- 
tems in the Galaxy as the function of the common 
envelope stage efficiency acE- 
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Fig. 3. — This figu re shows OCCO criterion 
(|Lipunov et al. Ill996bh for the ratio ^^"'+"^ 



Nf 



is the characteristic kick velocity. "Width" of the 
filled area depicts various values of acE in the 
range between 0.2 and 1.0. Observational value of 
the ratio ^f^f^ is ~ 0.001. 



500 600 700 800 



Fig. 2. — This figure shows how the ratio 
depends on the kick velocity vq. Here Np^r+NS is 
the calculated number of binary neutron stars with 
radio pulsars and Npsr is the calculated number 
of all radio pulsars, binary and single. "Width" of 
the filled area depicts various values of ac e in the 
range 0.2 - 1.0. 
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Fig. 4. — Birth frequency for different types of X-ray sources in the Galaxy. Marks in the Figure are: 1, 
NA+I; 2, NA+II; 3, NA+III; 4, NA+Be; 5, BH+II; 6, BH+III; 7a, SNA+III, collimation angle a = 10°; 
7b, SNA+III, a = 1°; 8a, SBH+III, a = 10°; 8b, SBH+III, a = 1°; 9a, SBH+III, a = 10°; 9b, a = 1°. For 
curves 9a, 9b minimal initial mass of the primary star is I2OM0, in other cases it is equal to lOM©. These 
calculations were conducted using stellar wind type A. 
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Fig. 5. — Birth frequency for different types of X-ray sources in the Galaxy. Marks in the Figure are: 1, 
NA+I; 2, NA+II; 3, NA+III; 4, NA+Be; 5, BH+II; 6, BH+III, 7a, SNA+III, coUimation angle a = 10°; 
7b, SNA+III, a = 1°; 8a, SBH+III, a = 10°; 8b, SBH+III, a = 1°; 9a, SBH+III, a = 10°; 9b, a = 1°. For 
curves 9a, 9b minimal initial mass of the primary star is I2OM0, in other cases it is equal to lOM©. These 
calculations were conducted using stellar wind type C. 
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log L ^, ergs 

Fig. 6. — Cumulative luminosity functions of different types of X-ray sources in the Galaxy. Marks in the 
Figure are: 1, NA+I; 2, NA+II; 3, NA+III; 4, NA+Be; 5, BH+II; 6, BH+III; 7a, SNA+III, collimation 
angle a = 10°; 7b, SNA+III, a = 1°; 8a, SBH+III, a = 10°; 8b, SBH+III, a = 1°; 9a, SBH+III, a = 10°; 
9b, a = 1°. For curves 9a, 9b minimal initial mass of the primary star is I2OM0, in other cases it is equal 
to IOMq. These calculations were conducted using stellar wind type A. 
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Fig. 7. — Cumulative luminosity functions of different types of X-ray sources in the Galaxy. Marks in the 
Figure are: 1, NA+I; 2, NA+II; 3, NA+III; 4, NA+Be; 5, BH+II; 6, BH+III, 7a, SNA+III, collimation 
angle a = 10°; 7b, SNA+III, a = 1°; 8a, SBH+III, a = 10°; 8b, SBH+III, a = 1°; 9a, SBH+III, a = 10°; 
9b, a = 1°. For curves 9a, 9b minimal initial mass of the primary star is I2OM0, in other cases it is equal 
to IOMq. These calculations were conducted using stellar wind type C. 
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Fig. 8. — Cumulative luminosity functions of all investigated systems in the galaxy like the Milky Way. See 
Table 1 for numerical data. In this Figure a is the collimation angle in supercritical regimes of accretion. 
These calculations were conducted using stellar wind type A. 
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Fig. 9. — Cumulative luminosity functions of all investigated systems in the galaxy like the Milky Way. See 
Table 1 for numerical data. In this Figure a is the collimation angle in supercritical regimes of accretion. 
These calculations were conducted using stellar wind type C. 
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Fig. 10. — Birth frequency (a) of all investigated systems (differential function) and differential luminosity 
function (b) of X-ray binary sources in the Galaxy. Stellar wind type A. Marks in the Figure are: 1, 
collimation angle (for super critical regimes) a = 10°;2, a = l°. 
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Fig. 11. — Birth frequency (a) of all investigated systems (differential function) and differential luminosity 
function (b) of X-ray binary sources in the Galaxy. Stellar wind type C. Marks in the Figure are: 1, 
collimation angle (for super critical regimes) a = 10°;2, a = l°. 
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Fig. 12. — Cumulative luminosity functions of all investigated systems in the "elliptical" galaxy after the 
star formation burst. See Table 1 for numerical estimations. The curves in the Figure represent the next 
models: 1, stellar wind type A, collimation angle (for super critical regimes) a = 10°; 2, wind A, a = 1°; 3, 
wind G, a = 10°; 4, wind C, a = 1°. The time ranges after the star formation burst in the Figure are: a, 
0-10 million years; b, 10-100 million years; c, 100 million - 1 billion years; d, 1-10 billion years. 
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Fig. 13. — Evolution of the X-ray luminosity after the star formation burst (T — 0) in the galaxy with mass 
IO^^Mq. See Tab le 2 for numerical data. I n this Figure: 1, our calculations, stellar wind type A; 2, the 
result obtained bv iTatarinzeva et al. I ( 1989^ : 3, our calculations, stellar wind type C. 
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Table 1: Numerical approximation of the cumulative luminosity function in the spiral galaxy. Stellar wind 
type A. See Figure [5] for graphical data. 



Luminosity range, k' 
logL^, erg s^i 



31.0 ~ 


32.5 


-0.25 


32.5 - 


35.6 


-0.1 


35.6 - 


37.2 


-0.25 


37.2 - 


38.0 


-0.7 


38.0 - 


38.3 




38.0 - 


38.5 


« -8^= 


38.3 - 


39.2 


« -0.05'' 


38.5 - 


40.2 


« -0.05^= 


39.2 - 


41.1 


-0.7'' 


40.2 - 


42.2 


-0.7== 



"fit curve is N{> L) ^ L'' . 

'coUimation angle (for super critical regimes) a = 10° . 
'^coUimation angle (for super critical regimes) a = 1°. 



Table 2: Numerical approximation of the cumulative luminosity function in the spiral galaxy. Stellar wind 
type C. See Figure [9] for graphical data. 



Luminosity range, 
log La;, erg s~^ 



31.0 


- 32.5 


-0.4 


32.5 


- 35.5 


-0.1 


35.5 


- 37.2 


-0.3 


37.2 


- 38.0 


-0.6 


38.0 


- 38.2 


-1.8'' 


38.0 


- 38.2 


-3.6^= 


38.2 


- 39.1 


-0.1'' 


38.2 


- 40.1 


-0.3^= 


39.1 


- 39.5 


-3.5'' 


40.1 


- 40.5 


-3.5^= 


39.5 


- 41.0 


-0.75'' 


40.5 


- 42.0 


-0.75^= 



"fit curve is Af(> L) ~ L'=. 

''coUimation angle (for super critical regimes) a = 10°. 
'^coUimation angle (for super critical regimes) = 1°. 
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Table 3: Numerical approximation of the X-ray luminosity of the galaxy after the star formation burst. 
Stellar wind type A. 



Time rang 


:e, 




pa 


10*5 yr 








4-10^-1 


- 103 


3 - 10^^ 


-2.5 


1-103-2 


- 103 


3.6 - 10^^ 


-1.56 


2 - 103 - 2, 


.5 • 103 


2 • 10^2 


^4 


2.5 - 103 - 


1 - 10^ 


3 - 10*^ 


-1.8 



"fit curve is L{T) = ci{T/10'^yr)P erg s"!. 



Table 4: Numerical approximation of the X-ray luminosity of the galaxy after the star formation burst. 
Stellar wind type C. 



Time range, 




pa 


10*5 yr 






4-10-1-40 


1.5 - 10*^ 


Ri 


1 - 40 - 1.5 - 103 


2 - 10*5 


-2.7 



"fit curve is L{T) = ci{T/10'^yr)P erg s"!. 
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